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Charging graphene nanoribbon quantum dots 
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We describe charging a quantum dot induced electrostatically within a semiconducting graphene 
nanoribbon by electrons or holes. The applied model is based on a tight-binding approach with 
the electron-electron interaction introduced by a mean field local spin density approximation. The 
numerical approach accounts for the charge of all the pz electrons and screening of external poten¬ 
tials by states near the charge neutrality point. Both a homogenous ribbon and a graphene flake 
embedded within the ribbon are discussed. The formation of transport gaps as functions of the 
external confinement potential (top gate potential) and the Fermi energy (back gate potential) are 
described in good qualitative agreement with the experimental data. For a hxed number of excess 
electrons we find that the excess charge added to the system is, - depending on the voltages defining 
the work point of the device: i) delocalized outside the quantum dot - in the transport gap due to 
the top gate potential ii) localized inside the quantum dot - in the transport gap due to the back 
gate potential or Hi) extended over both the quantum dot area and the ribbon connections - outside 
the transport gaps. The applicability of the frozen valence band approximation to describe charging 
the quantum dot by excess electrons is also discussed. 


I. INTRODUCTION 

The dispersion relation of graphen^ - gapless and 
linear near the charge neutrality point - p revents elec¬ 
trostatic confinement of charge carriers^. The energy 
gap in the band structure is opened for finite strips 
of graphen^SHZl which allows for carrier confinement by 
external potentials defining quantum dots inside the 
na norib bonAll nanoribbons exhibit a transport 
gapP^ near the neutrality point. Formation of the trans¬ 
port gap can be described theoretically already for non¬ 
interacting electrons as due to lateral confinement and 
edge disordeJ^. Nevertheless, the Anderson localization 
effects are reinforced by the electron-electron interaction 
and according to the present understanding, the trans¬ 
port gap is stabilized by a spontaneous formation of mul¬ 
tiple Coulomb islands along the ribbon. The current 
peaks that appear for a biased ribbon inside the trans¬ 
port gap are a ttributed t o electrons hopping between the 
charge puddlej^ ^ l ^^ l ^^ l ^^^ ^^. Charging of the intentionally 
define d qua ntum dots appears also within the transport 
gappHnsiSo]^ The spontaneous disorder-induced quantum 
dots and the ones defined electrostatically are charac¬ 
terized by Coul omb diamonds in the charge stability 
diagram j^^ l ^^ l ^^ l In part of the studies t he nanorib bons 
contain an additional inline graphene fiak^I^E^l^SIllHlll in 
the gated region, which is useful in separatioiP^ of the 
charging effects due to the intentionally introduced quan¬ 
tum dot and the spontaneous Coulomb islands along the 
ribbon. 

The electronic properties of quantum dots made of fi¬ 
nite graphene flakes have been extensively studied in a 
number of paperP^EH ^nd only recently a finite graphene 
flake was successfully connected directly to source-drain 
electrodes without nanoribbons feeding the charge to the 
fiak^^. Localization of electrons within a quantum dot 
induced electrostatically inside a graphene sheet with the 


mass gap induced by the substrate have been studied in 
Ref. [38l and the one-electron properties of these dots 
were discussed in Refs. 13^01 

The purpose of the present paper is a simulation of 
charging a quantum dot defined by an external poten¬ 
tial within a graphene nanoribbon. We consider both 
homogenous ribbons and the ones containing an inline 
flake of graphene. We prepared a numerical model 
that covers formation of both n-type and p-type quan¬ 
tum dots confining conduction band electrons or valence 
bands holes depending on the sign of the external poten¬ 
tial. For the purpose of the present study we consider 
a semiconducting armchair nanoribbon using the tight- 
binding approach for the single-electron spectrum and 
the LSDA potentials describing the electron-electron in¬ 
teraction. The model accounts for contribution of all 
the Pz electrons within the sample in formation of the 
potential profile along the system. We find that within 
the transport gap the quantum dot charging in terms of 
the chemical potentials can be very well described in a 
single-electron basis limited to an energy range near the 
neutrality point. The charge density of the states deep 
inside the valence band can then be considered frozen 
without any detectable deviation from the electron-hole 
symmetry in terms of electron and hole charging pro¬ 
cesses. We discuss the transport gap and the charge dis¬ 
tribution within the entire system including an integer 
quantization of the quantum-dot-confined charge. For a 
neutral system the modification of the lateral confine¬ 
ment - with introduction of the flake into the ribbon - 
does not induce formation of a localized charge puddle. 
Nevertheless, the excess electrons - when present within 
the system - get localized within the flake already in the 
absence of the external potential. A good qualitative cor¬ 
respondence with the transport gaps observed in Ref.l^^ 
as functions of the back gate and the top gate potentials 
is found. Applicability of the approximation of a frozen 
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valence band for charging the dot with excess electrons 
is also discussed. 
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FIG. 1: (a) Schematics of the ribbons considered in this work 
with and without the central flake. The considered nanorib¬ 
bons have armchair edges and the number of atoms across the 
channel (18 and 36) corresponds to a semiconducting disper¬ 
sion relation of pz tight-binding Hamiltonian. The length of 
the ribbon is = 61 nm. The hexagonal flake of width 
Lx 2 = 16.7 nm is embedded inside the ribbon the width 
Lxi = 2.15 nm (left plot) We consider also homogenous rib¬ 
bons of width Lxi = 4.3 nm (left plot), (b) External potential 
defining the quantum dot as given by Eq. ( |12| for W = 0.5 
eV. 


II. THEORY 

We consider a finite ribbon with armchair edges [see 
Fig. [^a)] with or without a hexagonal flake defined in 
its center. We use the single-electron Hamiltonian as 
given by the tight-binding method for Pz orbitals 

H tc\cp+h.c., ( 1 ) 

where t = —2.7 eV is the hopping parameter, and c] is 
the operator of electron creation at i-th atom and the 
summation runs over nearest neighbor atoms. 

Upon diagonalization of the Hamiltonian we obtain the 
eigenstates of form 


Na 

= ^CijP^{Y - ri). (2) 

i=l 

We assume 18 or 36 atoms across the ribbon, so that the 
system - in the Pz tight-binding approach has a semicon¬ 
ducting dispersion relation [see Fig. 1]. The charge re¬ 
distribution that we describe requires application of large 
external potentials, for which the potential at the ions is 
far from the charge neutrality point and the linear dis¬ 
persion relation of graphene. For this reason we need 
to employ an atomistic tight-binding approach instead 


of its continuous low-energy approximation in the Dirac 
form. The atomistic approach has its price on computa¬ 
tional complexity and memory consumption, that limits 
the size of the systems than can be studied. Neverthe¬ 
less, the obtained picture is qualitatively independent of 
the size the systems and the ones considered in this work 
remain within an experimental reach. For applications of 
the graphene nanoribbons in electronics their widt h be¬ 
low 10 nm is needed to switch off the currenfPIEll. The 
external potential well defining the size of the quantum 
dot is 15-20 nm [Fig. 1(b)], while the estimated diameter 
of the quantum dot in a gated system is 25-40 nnJ^. 

The system of interacting electrons is treated by a 
mean field approach. The number of ions in the con¬ 
sidered systems and the number of Pz orbitals in the ba¬ 
sis ranges up to about = 16 thousands, which sets 
the size of the basis (§. The considered systems are 
close to the charge neutrality, with the number of elec¬ 
trons within the ribbon close to the number of positive 
ions (net charge ±24e). We have found that the ground- 
state properties of the system can be quite accurately 
described using a basis of a few hundred wave functions 
of type for the energies near the neutrality point. 

We order the single-electron basis functions by the en¬ 
ergy eigenvalues Sj , where j is an integer non-zero index 
j G [—Na/2^Na/2] {Na is an even number). For the 
LSDA^^ calculations we use a basis of form 

m 

<t>kiy) = ( 3 ) 

j = -m 

where a is the spin index. The basis (§ contains m 
single-electron states of both the valence (negative j) and 
the conduction bands (positive j index). The states that 
are deep below the energy gap with j G [—Na/2^ —m — 1] 
are considered frozen. Their contribution to the electron 
spin density is calculated once and for all after diagonal¬ 
ization of the single-electron Hamiltonian 

—m—1 

= Y ( 4 ) 

j=-N^/2 

For the simulation of the mean-field potential we neglect 
the overlap of neighbor orbitals; the approach is equiva¬ 
lent to two-center treatment of the Coulomb integralJ^. 
Accordingly, the contribution of the j— th orbital to the 
spin density at /-th ion is given by \cijPz{r — rj)p, and 

—m—1 —m—1 

i=-JV„/2 j=-N^l2 

( 5 ) 

The electron-electron interaction is introduced by the 
LSDA mean field potentiaP^, with the Hamiltonian of 
form 


Hdft — 
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3d' 


(6) 
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where is the electron creation operator in the eigen¬ 
state j of Hamiltonian (1) with spin a. The DFT poten¬ 
tial enters the hopping parameters 

where n^(r) is the spin a density 

n^(r) =n'^(r) + n^i^(r). (8) 

In the above formula is the spin density of elec¬ 
trons occupying energy levels close to the neutrality point 
which are determined by DFT, 


where the summation runs over the ions i, and n(ri) = n^. 
We use the didectric constant e = 6 as for the graphene 
grown on 


In Eq. (11) Vc is the potential of the carbon ions which 


is calculated in a similar manner 


Vc{vi) = - 


47reeo 






rj| 




n'"(r)= ^ f{Ep,E^,)mr)(^ 


( 9 ) 


(14) 

The exchange and correlation potentials are taken ac¬ 
cording to the Perdew-Zunger parametrizatioiP^ for the 
spin density given at the ions. Once, the potential are 
known, we calculate the hopping elements Q 


k=—m 


Na 


where is the Fermi-Dirac distribution for the 

Fermi energy Ep^ and E^ are the eigenvalues of the DFT 
energy operator (H). The Fermi energy Ep is found from 
the normalization for the electron density n(r) = n^(r) -f 
n'^(r). 


= V Wc* 

JJ / ^ ^ K 




(15) 




/ 


where n^, 

We evaluate the chemical potential of the N electron 
system 


n(r)(i^r = 


( 10 ) 


where N is the number of electrons within the system 
N = Ni-\-Ne^ with Ni that stands for the number of ions 
within the nanoribbon, and for the number of excess 
electrons with respect to a charge neutral system. 

The potential defining the hopping parameters 0 
is given by 

W(r,n^(r),nHr)) = V.t(r)+n/(r, n)+l/c(r)+Wc(«t n^-). 

The external electrostatic potential in Eq. ( [ll| ) defining 
the quantum dot in the center of the ribbon is modelled 

by 

V.t(r) = -W^exp(^-(|)"^, (12) 

where 2R = 20 nm is the length of the quantum dot - 
see Eig. [^b) for the profile plotted for W = 0.5 eV. 
The electrostatic potential defining the quantum dot is 
assumed with an axial symmetry, since the gates are usu¬ 
ally applied perpendicular to the nanoribborP^. The po¬ 
tential as induced electrostatically is smooth and - due 
to the screening- short range. The exact profile of the 
confinement potential depends on the size of the gates 
and their distance to the area occupied by electronP®^. 
The qualitative features for the charging of the quantum 
dots remain independent of the potential profile, as long 
as it is smooth and short range. The results for a mod¬ 
ified form of the potential is given in Section IV.E, with 
power of 2 instead of 4 in Eq.(12). 

The Hartree potential l///(r,n(r)) = / , 

is evaluated as 


,, _ rptot rntot 

iiN — ~ 

where E]^^ is the total energy, calculated as 


= r + E, 

with the kinetic energy 


■ Ecr. 


Eh, 


(16) 

(17) 


T= -- 


2mo 


fit p 


cr=t,t k=—m 


(18) 

which can be expressed using the single-electron energy 
eigenvalues Sj for operator (1) with eigenfunctions ^ 


^ = E E E 1^. 


'jk I 


-J' 


(19) 


cr="|"4^ k= — m 


j = -m 


Eext,c in Eq. 0 is the contribution to the energy of 
the external potential and the carbon lattice 


Eext,c = J (Vextir) + Vc{r))n{r)d^ 


( 20 ) 


Exc in Eq. (17) is the Perdew-ZungeiS^ exchange- 
correlation energy and the last term of ( p!7| ) is the Hartree 
energy 


Eh = X 


= 11 VH{r)n{r)d^ 


( 21 ) 


In the summation of the single-electron contributions to 
the kinetic energy T (18) we account only for the 2m or- 


Vh{vi) = 


47reeo 


E 


rii 


■i-ii 


ni 


J 


|r-rz| 


(13) 


bitals forming the basis for the DET calculation, since the 
contributions of lower energy levels - considered frozen 
- cancel in the evaluation of the chemical potential any¬ 
way. However, in the exchange-correlation and Hartree 
potential, the total electron density, including the frozen 
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orbitals needs to be taken into account. For convergence 
of the DFT equations we use the BroydeiP^ method. The 
convergence is further additionally enhanced by anneal¬ 
ing of the temperature. We start from 20 K and take 
down the temperature to T = 10 K. This range of the 
temperature corresponds to the one actually used in the 
experimental studies of graphene quantum dotP^. In the 
model system, which is small, below 10 K no further po¬ 
tential dependence on the temperature is observed. 


III. RESULTS AND DISCUSSION 
A. Single-electron spectra 

Let us first characterize the single-electron spectrum 
in the absence of the electron-electron interaction for the 
quantum dots induced within the ribbon. The energy 
spectrum is displayed in Fig. [^as a function of the quan¬ 
tum dot potential for a nanoribbon of length L = 61 nm 
and 18 or 36 atoms across the ribbon - Fig. ia) and 
(b), respectively. The colors of the plotted energy lev¬ 
els indicate the localization of the corresponding states. 
The red energy levels are entirely localized within the 
quantum dot area [—1.2R, 1.2R]. For the quantum dot 
induced by the external potential in the homogenous rib¬ 
bon the localized energy levels appear for VF 7 ^ 0. The 
dot-localized (red) energy levels for kF > 0 correspond to 
localized states of the conduction band (n-type quantum 
dot) and for kF < 0 to the localized states of the valence 
band (p-type quantum dot). The energy levels plotted in 
blue in Fig. [^are localized entirely outside the quantum 
dot, generally near the ends of the ribbon, and they are 
insensitive to the dot potential. 

At the charge neutrality point {E = 0) [Fig. ia,b)] 
one observes an energy gap between the dot-localized 
(red) energy levels as a function of the dot potential 
kF. For = 0 we find localized conductance (valence) 
band states appearing for increasing kF > 0 (decreasing 
kF < 0). A gap is also observed as a function of the en¬ 
ergy for a fixed kF. In the experimental conditions the 
transport gap is observed as a function of both the dot 
potential (top g^e potential) and the Fermi energy (back 
gate potential/^. The latter corresponds to variation of 
the energy for a fixed kF value. 

Figure l^c) shows the spectra for the hexagonal flake 
embedded inside the ribbon with 18 atoms across the 
channel. For the flake the number of localized (red) en¬ 
ergy levels [Fig. [^c)] within the gap is much larger as 
compared to the nanoribbon quantum dot [cf. Fig. [^a)], 
and the localized energy levels appear already at kF = 0. 


B. Charging the nanoribbon quantum dot 

In Fig. |^a,c) we plotted the chemical potentials for the 
homogenous nanoribbon with 18 [Fig. [^a)] and 36 [Fig. 
[^c)] atoms across the channel. The considered number of 



FIG. 2: Single-electron spectra for a quantum dot of length 
2R = 20 nm defined within a semiconducting armchair 
nanoribbon of length L = 60.1 nm and 18 (a), or 36 (b) atoms 
across the channel. The colors of the lines indicate the extent 
of the electron localization within the dot defined as the prob¬ 
ability to find the electron within the segment [—1.2R, 1.2R] 
(see the color scale at right-hand side of the plots). In (c) 
the results for the hexagonal flake inside the nanoribbon with 
18 atoms are given. The black insets symbolically denote the 
studied system. 
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FIG. 3: (a) Chemical potentials for Ne — —24,... 24 excess electrons in a 18-atoms wide nanoribbon as a function of the QD 
potential W. (b) Charge contained within the quantum dot (segment [— 1.2i^, 1.2R]) defined within the nanoribbon 18- atoms 
wide as a function of W for Ne = —24,... 24 excess electrons. (c,d) - same as (a,b) for nanoribbon with 36 atoms across the 
channel. (e,f)- same as (c,d) but for the hexagonal flake with 44 atoms along the side of the hexagon. In (a-f) and elsewhere 
in the paper the length of the system is L = 60.1 nm. (g,h) - same as (g,h) but with longer ribbons with L increased to 90 nm. 
Insets in the figures shows an enlarged fragment from the graph. The orange lines in (a) mark the transport gap due to the 
Fermi energy (solid lines) and (b) the external potential (dotted lines) - see text. The thicker line in (a) shows the chemical 
potential of 8 electron confined in the dots /is - which is discussed in detail in the text. 
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Rd (nm) Rd (nm) 

FIG. 4: (a,b) Charge contained within a segment of 

[—1.2Rd,1.2Rd] for Ne = 1,2,... excess electrons for the 
quantum dot defined within the 36-atoms wide ribbon width 
for W = 0.5 eV (a) and W = 0.8 eV. For the energy spec¬ 
trum and the charge contained within the quantum dot see 
Fig. |^a,b) for W = 0.5 eV (a) and W = 0.8 eV (b). (c,d) 
- same as (a,b) only for a semiconducting 18-atoms ribbon 
containing a hexagonal flake across the channel. 


excess electrons varies between = —24 and = +24. 
For a system in contact with an external electron reser¬ 
voir the ribbon will be charged by electrons up to the 
electrochemical potential of the reservoir /i^, i.e. the rib¬ 
bon will contain N electrons for /tat+i > Mr > 

A necessary condition for the current flow at low bias is 
Mr = M 7 otherwise the charge transfer across the system 
is in the Coulomb blockade regime. The horizontal green 
line in the Fig. [^a,c) marks Mr = —803.92 meV, which is 
the value of the Fermi energy calculated for the neutral 
ribbon Ne = 0. This value - as we find - is independent of 
the size of the ribbon and corresponds to the center of the 
gap at the chemical potential spectrum for VF = 0, and 
can be treated as the charge neutrality point for the inter¬ 
acting system. The chemical energy gap - correspon ding 
to the gap in the transport experiments^EoEMlIlol _ 
the ribbon with 36 atoms across is increased from 268.5 
meV [Fig. to 304.7 meV [Fig. [^c)] by the electron- 
electron interaction. The chemical potential spectra - 
calculated with the restricted basis - are perfectly sym¬ 
metric with respect to the point inversion across the cen¬ 
ter of the energy gap. The electron-hole symmetry found 
in this result is quite remarkable, given that the frozen 
orbitals correspond entirely to the valence band. 

Note, that the present calculation deals only with the 
ground state of the electron system with Ne excess elec¬ 
trons within the dot. The gap in Fig. [^a) appears in 
the chemical potentials for varied electron numbers and 
not between the excited states of the single-electron or 
quasiparticle spectrunP^H^. The chemical potential gap 
determines the charging properties of the system and not 


- for instance - the optical oneP^^^ for a fixed confined 
charge. 

In the experimental conditions the center of the chem¬ 
ical potential gap can be set above or below the electro¬ 
chemical potential of the external reservoir depend¬ 
ing on the back gate potential. For below (above) 
this value, the ribbon outside the QD displays a p-type 
(n-type) conductivity. On the other hand the exter¬ 
nal potential (12) forms a confinement for conduction 
band electrons for VF > 0 or for valence band electrons 
{W < 0). For jj^r > —803.92 meV and VF < 0 we have 
n-type conductivity in the ribbon outside the quantum 
dot and confinement of the electrons of the valence band 
(p-type) within the quantum dot. This case is labeled by 
’NPN’ in Fig. [^a,c,e,g). Other cases of confinement are 
denoted in a similar manner. 


The electrons added to the nanoribbon fill the energy 
levels up to the electrochemical potential of the external 
reservoir but do not necessarily charge the quantum 
dot itself. The charge confined in the quantum dot Qd 
- for the nanoribbon with 18 and 36 atoms across the 
channel- as a function of W is plotted in Fig. |^b,d) for 
Ne = —24, —23,... 24. The value of Qd is found by in¬ 
tegration of the electron charge within the [—1.2i?, 1.2R] 
segment of the nanoribbon upon subtraction of the posi¬ 
tive ions present within this region. With the green (red) 
line we plotted the number of excess electrons confined 
inside the dot for = 0 and p^ = —803.92 meV. The 
number of confined electrons jumps in steps between the 
lines corresponding to a given number of electrons within 
the nanoribbon. For p^ = 0 - the value above the chem¬ 
ical potential gap - the entire nanoribbon is filled with 
electrons, over the potential of the dot - see the red line 
Fig. §h ,d) - and the charge within the dot takes nearly 
continuous and not quantized values. On the other hand 
for fir = —803.92 meV - the center of the gap for VF = 0, 
the number of electrons confined within the dot takes 
only integer values (green line in Fig. [^b,d)) - and cor¬ 
responds to charging the quantum dot with subsequent 
electrons. Note, that the entire spectrum of chemical po¬ 
tentials in Fig. |^a,c) is slanted, i.e. the horizontal green 
line in Fig. [^a,^ gets closer to the valence (conduction) 
band for VF > 0 (W < 0). As the number of electrons 
inside the nanoribbon grows, so does the potential of the 
background electron charge, hence the slope of the spec¬ 
trum, shifting up the center of the gap for VF >0. 

For W = 0.5 eV the chemical potentials of Ne G (2, 5) 
excess electrons appear within the gap of the 36-atoms 
wide ribbon of Fig. |^c). The chemical potentials within 
the gap acquire a stronger dependence on W and the 
charge within the dot for these electron numbers [inset in 
Fig. [^d)] becomes locally independent of W. For W = 
0.8 eV the chemical potentials of Ne G [6, 9] electrons 
appear inside the gap [Fig. id)], and the charge near 
the dot approaches a constant value before the end of 
the ribbon [Fig. gb)]. For TVe < 6 the quantum dot 
attracts the electrons from the entire nanoribbon and for 
Ne > 9 the excess electron charge is spread all along the 
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FIG. 5: Excess electron charge (a,d,g) and the total potential (b,e,h) for the nanoribbon of Fig. |^a,b) at W = 0.5 eV for 
Ne = 1,5 and 9 respectively. (c,f,i) present the cross sections of a-b, d-e, g-h, along the axis of the ribbon, respectively, 
compared the external potential of Eq. (12). 


system. In order to inspect the correspondence between 
the charge localized near the dot, the excess charge in 
the nanoribbon and the position of the chemical potential 
within the gap in Fig. |^c) we displayed in Fig. [^a,b) the 
charge within the [—Rd^Rd] segment of the nanoribbon 
for a varied number of excess electrons as a function of Rd 
for potentials W = 0.5 eV and W = O.S eV, respectively. 
For low Ne we find that the charge in the neighborhood 
of the dot Qd exceeds the number of excess electrons (see 
Fig. I^a) for A^e = I 5 • • • 4 and Fig. [^b) for Ne = .. .7), 

and a local maximum of Qd as a function of Rd is formed. 
The excess electron density for Ne = I and W = 0.5 eV is 
displayed in Fig. [^a), with the total potential plotted in 
Fig. ||:b). The external potential defining the quantum 
dot attracts also the electrons occupying the states of 
the valence band. The electron charge within the dot is 
accumulated at the expense of the neighborhood of the 
nanoribbon which acquires a positive charge - see Fig. 
[^c). For Ne = I the positive background is left along 


the entire ribbon up to its slanted ends (Fig. ga)). For 
W = 0.5 eV the local maximum of Qd{Rd) is no longer 
present for Ne = [Fig. [^b)]. For Ne = the excess 
charge is present only within the quantum dot [Fig. |^d)] . 
Note, that for Ne = the electron system spontaneously 
forms barriers in the total DFT potential separating the 
quantum dot from the rest of the system [Fig. The 

appearance of the barriers explains the saturated excess 
charge observed as a function of Rd in Fig. I^a). For 
Ne > b the value of Qd grows monotonically with Rd up 
to the ends of the nanoribbon [see Fig. [^b)], indicating 
that the excess electrons do not fit inside the quantum 
dot [see Fig. [^g) for Ne = 9] and are only confined due 
to a finite length of the nanoribbon. In this case the total 
potential is no longer constant outside the quantum dot 
[Fig. gh,!) for iVe = 9]. 

Figurega) shows the charge localized within the quan- 
turn dot for a fixed VF as a function of The values 
plotted by thick lines correspond to the chemical poten- 














































FIG. 6: The charge localized in the quantum dot Qd of a 36- 
atoms-wide nanoribbon (a) and the hexagonal flake embedded 
within 18-atoms-wide ribbon (b) for W = 0, 0.5 and 1 eV as a 
function of the electrochemical potential of the external elec¬ 
tron reservoir jir- The region of the transport gap is marked 
with the vertical lines. In (b) with the lines (dots) we marked 
the results obtained for nanoribbon length L = 60.1 nm (90 
nm). For the chemical potential spectra see Fig. |^a) and 
Fig. [^e,g), for the system without (with) the flake inside the 
nanoribbon, respectively. 


a) 


c) 


30 

20 

10 

0 

-10 

-20 

-30 


W-0.5 eV 






-6-3 0 3 6 
X (nm) 


30 

20 

10 

0 

-10 

-20 

-30 


W=0'.6 eV 




-6-3 0 3 6 
X (nm) 


30 

20 

10 

0 

-10 

-20 

-30 


W=1.6 eV 


I 

I 


e) 


-6-3 0 3 6 
X (nm) 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 


b) 


E 

0.7 

0 . 6 ^^ 

0.5 

0.4 — 
0.3 

0.2 c 
0.1 « 
0 a 
- 0.1 


d) 



30 

20 

10 

0 

-10 

-20 

-30 


W=0 eV 




-6-3 0 3 6 
X (nm) 


30 

20 

10 

0 

-10 

-20 

-30 


I 

■ 

■ 

m 

m 




-6-3 0 3 6 
X (nm) 


30 

20 

10 

0 

-10 

-20 

-30 


W=2 eV 

H 


M 

M 


I 


f) 


-6-3 0 3 6 
X (nm) 


E 

0.8 

0.7 

0.6 Td 
0.5 

0.4 >, 
0.3 ■(/) 
0.2 0 
0.1 "o 
0 ^ 
-0.1 


03 


1.8 

1.6 

1.4 

1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 


0 

■D 

0 

03 

0 


0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 


0 

■O 

0 

03 

0 


FIG. 7: The charge density of 8-th excess electron added to 
the system calculated as a difference of the total charge for 
Ne = S and Ne — 7 for the ribbon of Fig. 3(a). The /xs line 
is marked by the thick solid line in Fig. 3(a). 


tial in the gap of Fig. ia). Only within the gap the 
charge acquires quantized integer values. 

In the experimenfP^ the transport gap is observed for 
both the back gate potential and the top gate potential 
- see the black cross in Fig. 1(c) of Ref. |T3l The gap 
as a function of the back gate potential^ corresponds to 
the gap observed in Fig. 3(a,c,e,f) as a function of the 
chemical potential (or /i^). We marked this gap by the 
solid orange lines in Fig. 3(a). The transport gap re¬ 
lated to the top gate potential [the diagonal black region 
in Fig. 1(c) or Ref. [13] is also present in our results. 
The gap due to the top gate corresponds to an extension 
of the central gap region in the direction parallel to the 
slopes of chemical potentials for the states localized in¬ 
side the quantum dot [dotted orange lines in Fig. [^a)]. 
Let us focus on the chemical potential of 8 excess elec¬ 
trons [bold black line in Fig. 3(a)] as a function of the 
gate voltage and look at the location of the 8 th electron 


added to the system as a function of the external (top 
gate) voltage. For this purpose we subtracted the total 
charge densities for Ne = S and = 1 systems for a var¬ 
ied number of electrons with the results plotted in Fig. 

For W = —0.5 eV [Fig. [^a)] we are in the center of 
the gap marked by the dashed orange lines in Fig. 3(a). 
The 8 th electron when added to the system occupies the 
ribbon connections and not the quantum dot. The flow 
of the current across the dot is then blocked, hence the 
vanishing current found in Ref. [13] for the transport gap 
defined by the top-gate potential. Note, that /ig within 
this gap have a negligible dependence onW - which also 
indicates that the 8 th electron is localized outside the 
quantum dot. When we leave the transport gap [Fig. 
[^b,c) for VF = 0 and W = 0.6 eV] the 8 -th electron oc¬ 
cupies both the ribbon connections and the quantum dot, 
allowing for the flow of the current across the system. In 
this region the chemical potential acquires a pronounced 
dependence on W. For W = 1 eV [Fig. 0 :d)] we are 
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in the center of the gap on the Fermi energy scale (back 
gate of Ref. [13]) and the 8-th electron occupies the dot 
only and the variation of /ig with the external potential 
is the strongest. In the region marked by the orange solid 
lines in Fig. 3(a) the transport across t he qu antum dot 
can occur via the dot localized states^lMl^, naturally 
provided that the ribbon connections have non-zero den¬ 
sity of states in this energy range. When we leave the 
transport gap [Fig. [^e,f) for W = 1.6 eV and 2 eV] the 
electron occupies both the dot and the ribbons allowing 
for the electron flow across the system. 

The transport gap due to the top gate potential in Fig. 
1 (c) of Ref.l^ depends strongly on the back gate poten¬ 
tial. An increase of the top gate potential in Fig. 1(c) of 
Ref.^ corresponds to an increase of W. A more negative 
back gate potential in Fig. 1(c) of Ref.^ corresponds to 
a more negative values of the chemical potentials in Fig. 
[^a). Hence for large positive W the gap appears for 
large negative chemical potentials in Fig. I^a) in a cor¬ 
rect qualitative correspondence with Fig. Re) of Ref.l^^, 
where the transport gap for large positive top gate po¬ 
tential appears for large negative back gate potentials. 


C. Charging the QD induced within the graphene 

flake 

Figure compares the chemical potential spectra 
and the charge localized within the quantum dot for a 
nanoribbon containing 18 atoms across the channel with 
[Fig. de ,f)] and without [Fig. [^a,b)] a hexagonal flake 
defined in the center of the ribbon. For the quantum dot 
induced within the flake the chemical potential gap as a 
function of W is much thinner as in the results obtained 
for non-interacting electrons cf. Fig. [ 2 ]; a) and Fig. [2];c). 
The dependence of the charge confined in the quantum 
dot as a function of W is displayed in Fig. [^f). For 
a fixed value of W the flake contains a larger number 
of electrons than for the ribbon QD [cf. Fig. [^f) and 
Fig. I^b)]. Nevertheless, no excess charge appears in the 
flake for a neutral system = 0 at W = 0 [Fig. [^f)]. 
In contrast to the nanoribbon QD - for Ne > 0 a few 
extra electrons are localized within the flake QD also at 
W = 0. In consequence, for the flake the central band of 
a quantized confined charge in Fig. if) is more or less 
constant as a function of W and in contrast to the 
nanoribbon [Fig. [^d)] where for W = 0^ Qd is only inte¬ 
ger for A^g = 0. The plots of the confined excess charge 
as a function of the cutoff-radius Rd for the flake quan¬ 
tum dot given in Fig. [^indicate a large flexibility of the 
confinement potential in terms of the number of confined 
electrons. For W = 0.8 eV [Fig. [^f)] a saturation of the 
excess charge is obtained for G (3,14) near Rd — 10 
nm. 

In order to estimate the effect of the finite size of the 
ribbons in the system we plotted in Fig. [^g) and [^h) 
the results for the hexagonal flake embedded in a longer 
ribbon. The total length of the system was increased 


from L = 60.1 nm - as in Fig. [^e,f) and elsewhere in 
this work, to 90 nm. We observe that the chemical poten¬ 
tials within the energy gap and the width of the energy 
gap are unchanged. Nevertheless, the chemical poten¬ 
tials which are outside of the gap are shifted towards 
the neutrality point - the band of chemical potentials 
for Ne G [—24, +24] electrons becomes thinner on the 
vertical scale. The longer ribbons reduce the interaction 
energy for the carriers within the ribbons, and the chem¬ 
ical potential approach the continuum limit. Concerning 
the charge confined within the dot displayed in Fig. |^h) 
- as compared to Fig. if) we can see that the central 
part of the figure - corresponding to a quantized charge 
within the dot - is left unchanged. Outside the energy 
gap the entire nanoribbon contains a larger number of 
electrons for a given /i^ and the line for /i^ = 0 - as 
compared to Fig. if) - tends to acquire a continuous 
dependence on W. Same conclusions are reached for the 
charge localized within the flake quantum dot for a fixed 
VF as a function of /i^ [cf. the solid and dashed lines in 
Fig. [^b) corresponding to L = 90 nm and L = 60 nm, 
respectively]. 


D. Frozen valence band charge 

The results presented so far were obtained within the 
method in which both the conduction and valence band 
states formed the basis for solution of LSD A equations. 
Let us suppose that we are only interested in the conduc¬ 
tion band states, or charging the quantum dot with elec¬ 
trons A'e > 0 and not holes - which appears for kF >0. 
Let us discuss the applicability of the method with the 
basis limited to conduction band states - for which the 
complexity of the calculations can be significantly re¬ 
duced. For this purpose we i) diagonalize Hamiltonian 
(1), fix the charge density (4) as obtained for a charge 
neutral system taking the summation up to j = —1 in Eq. 
(4), then Hi) the LSDA equations are solved in the basis 
given by Eq. (3) but with summation starting from j = 1. 
Since the Hamiltonian (1) does not contain any potential, 
the frozen charge density (4) of all the occupied valence 
band states below the energy gap is uniform and the same 
on each atom of the crystal lattice. Then, the solution of 
the DET equations is performed in the conduction band 
only, with the frozen valence charge producing only a 
constant potential shift, due to the exchange-correlation 
component of the Perdew-Zunger potential. 

The results for the chemical potentials for a nanorib¬ 
bon with 18 atoms across the channel are displayed in 
Eig. I^a). The basis limited to the conduction band cor¬ 
rectly reproduces the chemical potentials above and in¬ 
side the gap. No additional shift was necessary to make 
this data coincide. As W is increased, the chemical po¬ 
tential spectrum for low > 0 continues into the quasi¬ 
continuum of the chemical potentials of the valence band. 
The basis limited to the conduction band state predicts 
that the charge confined within the quantum dot remains 
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FIG. 8: (a) Chemical potentials for Ne = —24,... 24 excess 
electrons for the ribbon with 18 atoms across the channel as 
calculated by the standard method of the present work (red 
lines for Ne — 1,... 24 and grey lines for Ne — —24, ••• — !) 
and with the basis of conduction band only (black lines), with 
all the single-electron states below the energy gap considered 
as frozen, (b) Charge conhned within the quantum dot for 
both approaches as a function of W. 


quantized at integer values for large values of W [black 
lines in Fig. I^b)], while in the exact calculation [red lines 
in Fig. [^b)]the system enters the energy continuum for 
the p-type conductivity. Near the bottom of the gap we 
find that for a fixed value of N^ the chemical potential 
calculated with the limited basis goes down on the energy 
scale as compared to the exact calculation [Fig. |^a)]. In 
this region we find that for a fixed Ne the exact calcu¬ 
lation gives values of the quantum-dot-confined charge 
which are i) larger [Fig. [^b)] than the ones obtained 
with the limited basis and ii) continuous as a function 
on W in contrast to the quantized values for the limited 
basis. This discrepancy is a result of filling the quantum 
dot by the electrons coming from the valence band from 
outside the quantum dot, which is present in the exact 
calculation but neglected in the limited basis. In the 
limited basis the excess charge confined within the quan¬ 
tum dot is underestimated, hence the underestimate of 
the interaction energy and in consequence lower values of 
chemical potentials as compared to the exact calculation 
[see the red and black lines in Fig. [^a) before the black 
lines enter the gray band of the valence band states]. 


E. Potential profile and edge disorder 

In this subsection we present the results for a mod¬ 
ified external potential profile and the disordered edge 
within the flake. The choice of the form of the external 
potential defining the quantum dot [Eq.(12)] was taken 
arbitrarily, as a smooth short-range potential. In order 
to demonstrate that the adopted results remain quali¬ 
tatively unchanged when the potential profile is varied 
we performed calculations for a Gaussian confinement 

= —W exp i-e. for a power of two in¬ 

stead of 4 in the exponent. The results for the chemical 
potential are displayed in Fig. to be compared with 
Fig. 3(a) for the potential given by Eq. (12). The pat¬ 
tern of chemical potentials remains unchanged, and only 
quantitative differences are observed. 

Eor the edge disorder effects we took the hexagonal 
flake with pure armchair edges end removed two consec¬ 
utive atoms every 10 of atoms along the edge. The form 
of the boundary is displayed in Eig. [lQ][a). The edge 
contains fragments zigzag edge and single separated ions 
within armchair fragments. The results with ideal arm¬ 
chair edges of the flake were presented in Eig. 3(b,c) for 
the chemical potentials and the charge localized within 
the dot. Again only quantitative differences can be spot¬ 
ted with no overall change to the quantitative description 
of the effect. 



FIG. 9: Ghemical potentials for Ne = —24,... 24 excess elec¬ 
trons for the ribbon with 18 atoms across the channel for 
conhnement potential —W exp ( — (^) ) instead of Eq. (12). 

For the results with Eq. (12) see Fig. 3(a). 


IV. SUMMARY AND CONCLUSIONS 

We have studied the charge distribution in a semicon¬ 
ducting graphene nanoribbon with a quantum dot de¬ 
fined in its center by an external potential. The calcu¬ 
lations were performed within the tight-binding method 
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with the electron-electron interaction introduced by the 
mean field DFT approach. 

We determined the excess charge localized within the 
dot as a function of the Fermi energy and the external 
potential within the transport gap and outside, when 
the entire structure is flooded by carriers. The calcu¬ 
lated chemical potential spectra bear signatures of the 
transport gaps which agree with the ones observed ex¬ 
perimentally as functions of both the back and top gate 
potentials. The inline flake embedded in the gated re¬ 
gion of the nanoribbon allows for confinement of the ex¬ 
cess charge already in the absence of the external poten¬ 
tial, however no excess charge within the flake is found 
for a charge-neutral system unless an external potential 
is present. We found that for the frozen valence band 
charge the charging of the quantum dot with excess elec¬ 
trons A^e > 0 can be quite well described for A^e > 0 
above the lower limit of the transport gap. However, the 
lower limit of the transport gap is overlooked by the basis 
limited to the conductance band states. Moreover, the 
limited basis neglects the QD charging with the valence 
band electrons which appears for large external poten¬ 
tials. The neglect leads to an underestimate of both the 
confined charge and the chemical potential, particularly 
near the lower limit of the transport gap and below on 
the Fermi energy / potential scale. 
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